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Abstract: We evaluate for arbitrary even dimensions the classical continuum limit 
of the lattice axial anomaly defined by the overlap-Dirac operator. Our calculational 
scheme is simple and systematic. In particular, a powerful topological argument is 
utilized to determine the value of a lattice integral involved in the calculation. When 
the Dirac operator is free of species doubling, the classical continuum limit of the 
axial anomaly in various dimensions is combined into a form of the Chern character, 
as expected. 
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In this note, we evaluate for arbitrary even dimensions d — 2n the classical con- 
tinuum limit of the lattice axial anomaly defined by the overlap-Dirac operator [1]. 
This quantity has been studied many times [2]-[8]; but all treat d — 4 case. There 
also exist general arguments [3, 9, 7, 8] that the classical continuum limit reproduces 
the axial anomaly in the continuum theory, when the Dirac operator is free of species 
doubling. If it is possible, however, it is certainly preferable to demonstrate this by 
an explicit calculation. We do this for arbitrary dimensions and obtain the expected 
result; the classical continuum limit of the axial anomaly in various dimensions is 
combined into a form of the Chern character [10] [see eq. (41)]. Our calculation 
consists of two steps: First we determine 0(a°) terms in the axial anomaly with- 
out expanding it with respect to the gauge coupling; this scheme is similar to that 
of refs. [11]— [14], [4, 5, 15]. Next, we evaluate a lattice integral involved; here a power- 
ful topological argument is utilized to determine its value (closely related topological 
arguments can be found in refs. [16, 17, 2]). A combination of these two schemes pro- 
vides a quite simple and systematic method to evaluate the classical continuum limit 
of the axial anomaly. To illustrate this point, we present a calculation of the axial 
anomaly with the Wilson-Dirac operator for arbitrary dimensions in appendix A. 

The overlap-Dirac operator [1] is defined by 1 

D = - \l - A{rf A)' 1 ' 2 ] , A = m -aD w , (1) 

A, = \ + V„) - arV^V J , (2) 

where m and r are free parameters. We assume r > in what follows. In the absence 
of the gauge field, the Dirac operator is free of species doubling if < m /r < 2. V M 
and V* in this expression are forward and backward covariant difference operators 
respectively 2 

V^(x) = -[U (x, Mix + afi) - rj}{x)\ , (3) 

CL 

V*iJ)(x) = - 4>(x) — U{x — ap,, /i)fy(x — ajl) . (4) 

(X 

By construction, the overlap-Dirac operator satisfies the Ginsparg- Wilson relation 
^ d+ iD + D^ d+ i = aD^d+iD [18]. Due to this relation, the lattice fermion action Sf = 
a d J2 x ' l P{ x )D'4 , (x) is invariant under the modified chiral transformation [19] 

6iJ>(x) = e-f d+1 (l - ^aD^j ip(x), Sip(x) = (l - ^ aD ^j Jd+it- (5) 



1 Out notations: a denotes the lattice spacing. Greek letters, /x, v, . . . run from 1 to d = 2n. 
Repeated indices are understood to be summed over, unless noted otherwise. {7^,7^} = 26^, 

it = 7p and Jd+i = (-«) n 7i ■■•Id; 7d+i = 1 and = follow from this. 

2 The link variable U(x, /j) denotes a matrix in a (in general reducible) unitary representation of 
the gauge group, to which the fermion belongs. 
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The fermion integration measure, however, acquires a non-trivial jacobian under the 
transformation [19] 



6Y[dip(x)dip(x) = -2 



Y\_ dip (x) dip (x), 



where 



q(x) = -jtr 



ld+i ( 1 - -aD(x,x) 



(6) 



(7) 



This framework thus realizes a desired breaking pattern of the chiral Ward-Takahashi 
identity even with finite lattice spacings; 3 it is especially suitable for a study of 
phenomena associated to the axial anomaly, such as the U(l) problem [20]-[23]. 
The lattice axial anomaly q(x) may also be regarded as the index density, because 
a d ^2 x q(x) 1S an integer depending on the gauge-field configuration [24]. 

We first note that a diagonal element of the kernel of an operator on the lattice 
can be expressed as 



O 



/7r r\du 



* d d k 
-* (27r) d 



e -~ikx/aQ e ikx/a 



(8) 



So we note 



{ ^ } e^ a f(x) = e ik *' a |t7m(^m + oQ„) + m + r £(c M - 1) + aR 



where 



and 



Sfj, — sin hp, — cos k^, 



From eqs. (1) and (7), we thus have 



1 



m, (9) 

(10) 

(11) 



q(x) 



d d k 



tr 



ld+i { -l^is^ + aQ^) +m + r 



2a d J-n {27i) d 

x —{is v + aQ v ) 2 + < m + r 



£(c„ ~l) + aR 



t 1 

in2 



J2(cu -l) + aR 



-1/2. 



-[lv,lp][Qv,Qp] + a 2 r[n/ u Q u ,R] 



1, (12) 



3 Note that the jacobian is unity for flavor non-singlet chiral transformations for which tre = 0. 
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where we have used tr 7^+1 = and r ) u r ) p = S up + [y v , 7 p ]/2. In this form, an expansion 
of q(x) with respect to the lattice spacing a is straightforward. In the classical 
continuum limit, the gauge potential A^(x) introduced by 



U (x, fx) —V exp 



a / dt AJx + (1 — t)ap) 
Jo 



(13) 



is regarded as a smooth function of x of the order 0(a°). The operators and R 
are of 0(a°) in this limit: 



and 



Qti = c M (9 M + Ap) + 0(a), no sum over //, 
R = Y l is li {d li + A li ) + 0(a), 



[Qv,, Qu] = c^CyF^ + 0(a), no sum over jj, and v, 
i?] = ic M ^ s^F Ml/ + 0(a), no sum over //, 



(14) 
(15) 



(16) 
(17) 



where F^ v = d^A u — d u A^ + [A^, A v \. Thus it is straightforward to find 0(a°) terms 
in eq. (12), because the trace over Dirac indices requires at least d gamma matrices: 



Namely, we have 

1 /■* d d A; /-1/2 N 



(18) 



/■7T 



(2tf) 



m o + ^^(c P - 1) 



2 s -n-1/2 



m + r 5Z( C cr - 1) 



+ tr<^7 d+ i7 a s 



n-l 



+permutations 



+ 0(a). (19) 



A little calculation using the identity (18) shows 
(-i) n /-1/2 



q(x) 



2(27r) fl! 

where the lattice integral I(m ,r) is given by 

I(mo,r) = / d d kl(k;m ,r), 



(20) 



(21) 
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with 



-n-l/2 



X(k;m ,r) 



II + 



m + ^^(Ci, - 1) 



X 



m o + ^ XX c /> - 1) + r E — 



and 
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A; M e M d 



7T , 37T 
< fc < — 

2 2 



(22) 



(23) 



Here we have shifted the integration region as A;^ e [— vr,7r] — > A; M e [— 7r/2, 37t/2] for 
later convenience. 

We now turn to a calculation of the lattice integral J(m ,r). We utilize the 
following topological argument. We first introduce a mapping from the Brillouin 
zone B to the unit sphere S d . The mapping is defined by 



O = + ^^(c M - 1), 
p 

6^ = Sp, for 



H = l,...,d, 



and 



0A 



(24) 
(25) 

(26) 



where xa (A — 0, 1, . . ., d) is the coordinate of in which the unit sphere Y^a x \ — 
1 is embedded. Since this mapping is periodic on the Brillouin zone B, we can regard 
B as a torus T d . Namely k^ — > defines a mapping / : T d — > S^. The crucial 
observation is that the volume form on this sphere coincides with the integrand 
of I (mo, r): 

1 



^ = ~d\ tA ^- A d XA ^ x M A • • • A drr A(j 

5^1 e Ao -A d 0A Q d9 Al A • • • A d6U d 



d! e° 

I(k] m , r) dki A • • ■ A dk d . 



(27) 



This shows that the integral of on a (sufficiently small) coordinate patch U on S d 
is given by 



Q — sgn 

Ju 



l(k j ;m ,r)} [ l(k; m , r) dfci A • • • A dk d , (28) 

where W (j — 1, . . . , m) is a component of the inverse image of C/ under /, / _1 (£7): 

f-\U) = U 1 U ••• Uf/ m C T d , (29) 
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and W (j = 1, . . . , m) is a certain point on [P. For a sufficiently small U, Uj are 
pairwise disjoint. We take preimages of a point y <E U under /, as fo 7 '. Then 

by summing both sides of eq. (28) over j, we have 4 

W l(k;m ,r)dh A---Adk d = (deg/) / ft, (30) 

~ Jltf ./[/ 

where the degree of the mapping / is given by 

deg/ = J2 sgn [j(F;m ,r)] . (31) 
/(fc J )=?/ 

In general, the degree of the mapping / : T d — > S d is defined by a sum of the signature 
of jacobian of the coordinate transformation between T d and S d over preimages of a 
point y G S^. An important mathematical fact [25] is that the degree takes the same 
value for all (regular) points y on S d . Thus, by summing eq. (30) over all coordinate 
patches U of S d , we have 

I(m ,r)= [ J( k; m ,r) dh A ■■■Adk d = (deg/) / ft, (32) 
where J sd ft is given by the volume of the unit sphere S d : 

[ n = vol(S d ) = 2d+1 f n] . (33) 
Js d d\ 

We may choose any point y on S d to evaluate the degree (31). We choose 
y — (1, 0, • • • , 0). This requires 

^(F) = ^ = 0, for /x = 1, . . . (34) 

and 

*o(F) = m ° + rE ^" 1) = 1. (35) 



m + rE (i (c fl - 1) 

Note that eq. (35) is equivalent to the condition m /r + X^(c M — 1) > 0. Now eq. (34) 
implies that k^ = or n for each direction /i. We denote the number of 7r's appearing 
in y by an integer n n > 0, 

^ = (vr,...,7r,0,...,0), (36) 



4 In deriving this relation, we have assumed that U is within the range of /, i.e., the inverse 
image / _1 (C7) is not empty. This relation itself, however, is meaningful even if U is not within the 
range of /, if one sets deg / = for such case. As a consequence, eq. (32) holds even if / :T d — > S d 
is not a surjection, i.e., not an onto-mapping. 
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irrespective of the position of 7r's. For a given n n , the number of such W is y n j. The 
second relation (35) on the other hand requires n n < m /(2r). At those k\ we have 



sgn 

Thus eq. (31) gives 



X(F;m ,r)l =I1 c m= (-!)""• (37) 



["io /2r] / j\ 

degf= £ (-!)"-■ (38) 

Combining eqs. (20), (32), (33) and (38), we finally obtain 

■ n [m /2r] / , \ 

= T^Tj E (- 1 )" [ n J tr • • • F^J + 0(a), (39) 

where we have used 

M/2 Mf41VD". (40) 



In particular, when < m /r < 2 with which the (free) overlap-Dirac operator is 
free of species doubling [1], eq. (39) reproduces the expected result. Multiplying the 
volume form d d x = dxi A • • • A dx d to q(x), we have 



q(x) d d x = trexp F 



+ 0(a), (41) 

d d x 



where the field strength 2-form F is defined by F = F^dx^ A dx u /2. This is our 
main result. When combined with cohomological arguments, it provides an absolute 
normalization of the lattice axial anomaly even with finite lattice spacings [26, 27] 
and with finite volumes [28]. 



A. Axial anomaly with the Wilson-Dirac operator 

With the Wilson-Dirac operator, the chiral symmetry is explicitly broken by the Wil- 
son term and the axial anomaly is regarded as arising from this explicit breaking [29]. 
The axial anomaly thus has a structure of a fermion one-loop diagram containing 
the Wilson term and, in our notation, it reads 



qw(x) = -rtr 



^ tr 



7, + iyrV;V,(AtA)- 1 At 



(A.I) 



where the parameter m in A of eq. (1) is replaced by am with m being the bare mass 
of the fermion. An expansion with respect to the lattice spacing is almost identical 
to the case of the overlap-Dirac operator in the text and we have 



q w {x) = 



(2tt)< 



Iw(r) e 



/Lill/l ■■■jll; 



(A.2) 
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where I w (r) = f B d d kI w (k;r) and 



1 2 



i) 



-n-l 



x < r 



1 2 



B 

L p 



1) 



+ ^ 2 E(cp-i)E- . (A.3) 



which is independent of the fermion mass. To evaluate this lattice integral, we again 
utilize the mapping (24)-(26) with mo = 0. The mapping — > x M (we do noi 
include Xq here) defines a mapping g : B — > where _B d is the d dimensional ball, 
£ d = {x M G R d | < !}• More precisely, we have to regard all points of the 

boundary of the ball, — 1, are identified, because fc = (0, ... ,0) is mapped 

to the boundary of the ball. The crucial relation this time is that the volume form 
of the ball, dx\ A ■ ■ • A drr^, coincides with the integrand of I w (r) as one can verify. 
These observations show 



Ur) 



(deg g) dxi A • • • A dx d = (deg g) vol{B d ) 

JB d 



(A.4) 



where the volume of the d dimensional ball, vo\(B d ), is given by vo\(B d ) = n n /n\ 
and the degree is 

deg#= sgn [j w (F;r)j . (A.5) 

g(ki)=y 

As the point y on B d , we choose the origin y = (0, ...,0). All momenta whose 
components k^ are either or n, except k = (0, . . . , 0), are mapped to y — (0, . . . , 0). 
Thus the degree is given by 



(A.6) 



Combining eqs. (A. 2), (A.4) and (A.6), we have q w (x) d d x = tr[expiF/(27r)]| d d :E + 
0(a) for arbitrary Wilson parameter r ^ 0. 
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